Effective free energy method for nematic liquid crystals in contact with structured 

substrates 
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We study the phase behavior of a nematic liquid crystal confined between a flat substrate with 
strong anchoring and a patterned substrate whose structure and local anchoring strength we vary. 
By first evaluating an effective surface free energy function characterizing the patterned substrate 
we derive an expression for the effective free energy of the confined nematic liquid crystal. Then we 
determine phase diagrams involving a homogeneous state in which the nematic director is almost 
uniform and a hybrid aligned nematic state in which the orientation of the director varies through 
the cell. Direct minimization of the free energy functional were performed in order to test the 
predictions of the effective free energy method. We find remarkably good agreement between the 
phase boundaries calculated from the two approaches. In addition the effective energy method 
allows one to determine the energy barriers between two states in a bistable nematic device. 

PACS numbers: 42.70.Df, 42.79.Kr, 61.30.Hn 



I. INTRODUCTION 

The interest in anchoring phenomena and phenomena 
in confined nematic liquid crystals has largely been driven 
by their potential use in liquid crystal display devices. 
The twisted nematic liquid crystal cell serves as an ex- 
ample. It consists of a nematic liquid crystal confined 
between two parallel walls, both providing homogeneous 
planar anchoring but with mutually perpendicular easy 
directions. In this case the orientation of the nematic 
director is tuned by the application of an external elec- 
tric or magnetic field. A precise control of the surface 
alignment extending over large areas is decisive for the 
functioning of such devices. 

Most studies have focused on nematic liquid crystals in 
contact with laterally uniform substrates. On the other 
hand substrate inhomogeneities arise rather naturally as 
a result of surface treatments such as rubbing. Thus the 
nematic texture near the surface is in fact non-uniform. 
This non-uniformity, however, is smeared out beyond a 
decay length proportional to the periodicity of the sur- 
face pattern. Very often the thickness of the non-uniform 
surface layer is considerably smaller than both the wave- 
length of visible light and the thickness of the nematic 
cell, i.e., the distance between the two confining paral- 
lel walls. Hence optical properties of the nematic liquid 
crystal confined between such substrates correspond to 
those resulting from effective, uniform substrates. 

More recent developments have demonstrated that sur- 
faces patterned with a large periodicity of some microm- 
eters are of considerable interest from a technological 
point of view (see, e.g., Ref. [l[ and references therein). 
A new generation of electro-optical devices relies on ne- 
matic liquid crystals with patterned orientation of the 
nematic director over small areas which can be achieved 
by chemically patterning the confining surfaces. For ex- 



ample, to produce flat-panel displays with wide viewing 
angles one can use pixels that are divided into sub-pixels, 
where each sub-pixel is defined by a different orientation 
of the nematic director, which is induced by the surface 
structure and subsequently tuned by the electric field. 
In addition to the technological relevance, nematic liquid 
crystals in contact with non-uniform substrates provide 
the opportunity to study basic phenomena such as ef- 
fective elastic forces between the substrates and phase 
transitions between various competing nematic textures 
(see, e.g., Ref. @ and references therein). 

Whereas the influence of homogeneous confining sub- 
strates on nematic liquid crystals is now well understood, 
the phase behavior of nematic liquid crystals in contact 
with chemically or geometrically patterned substrates is 
still debated. One might suppose that theoretical calcu- 
lations based on continuum theories should resolve the 
properties of nematic liquid crystals in contact with pat 
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However, such calculations are numeri- 
cally demanding because two- or three-dimensional grids 
have to be used because of the broken symmetry due 
to the surface pattern. Moreover, it is very challenging 
to determine metastable states and energy barriers be- 
tween them which are important for the understan ding 
of bistable nematic devices 133. [35 L |3 a . l37 l [38l l39l 140 l4ll . 
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In the present paper we adopt a different strategy which 
takes the advantage of the finite decay length character- 
izing the influence of the surface pattern on the nematic 
liquid crystal in the direction perpendicular to the sub- 
strate. We determine first an anchoring energy function 
and an average surface director orientation of the pat- 
terned substrate and obtain an effective free energy for 
the nematic liquid crystal cell under consideration. We 
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find remarkably good agreement between the phase dia- 
grams of various systems calculated using this effective 
free energy function on the one hand and the original free 
energy functional on the other hand. 



II. EFFECTIVE FREE ENERGY FUNCTION 

A. Continuum theory 

The continuum theory for liquid crystals has its ori- 
gin dating back to at least the work of Oseen [58[ and 
Zocher [59(. This early version of the continuum the- 
ory for nematic liquid crystals played an important role 
for the further development of the static theory and its 
more direct formulation by Frank 60] . The Frank theory 
is formulated in terms of the so-called nematic director 
n = n(r), |n| = 1, and its possible spatial distortions. 
The nematic director describes the direction of the lo- 
cally averaged molecular alignment in liquid crystals. In 
a nematic liquid crystal the centers of mass of the liq- 
uid crystal molecules do not exhibit long-ranged order. 
The molecules can translate freely while being aligned, 
on average, parallel to one another and to the nematic 
director. It is known that if an initially uniform nematic 
liquid crystal is distorted by external forces, it relaxes 
back to the uniform state after the disturbing influence 
is switched off, signaling that the uniform configuration 
represents a thermodynamically stable state. Therefore 
it is assumed that there is a cost in free energy associ- 
ated with elastic distortions of the nematic director of 
the form 



F e , as [n(r)] = \j v <Pr 



K n (V • n) 2 



+ K 22 {n • (V x n)) 2 + K 33 (h x (V x fi)) 2 ] , 

(1) 

where V is the volume accessible to the nematic liquid 
crystal and K\\, K22, and K33 are elastic constants as- 
sociated with splay, twist, and bend distortions, respec- 
tively. The elastic constants depend on temperature and 
are commonly of the order 10~ 12 to 1CP 11 N. Sometimes, 
for example, when the relative values of the elastic con- 
stants are unknown or when the resulting Euler-Lagrange 
equations are complicated, the one-constant approxima- 
tion K — Kn = K22 = K33 is made. In this case the 
elastic free energy functional reduces to 



F e i as [h(r) 



r (Vn) 2 



(2) 



In the presence of surfaces the bulk free energy F\, = F e i as 
must be supplemented by the surface free energy F s such 
that the total free energy is given by F = Ft + F s . In 
the corresponding equilibrium Euler-Lagrange equations 
5F/5h — 0, F s leads to appropriate boundary conditions. 
The description of the nematic director alignment at the 
surfaces forming the boundaries is called anchoring. In 



addition to the so-called free boundary condition where 
there is no anchoring, one considers weak and strong an- 
choring. If there are no anchoring conditions imposed on 
n at the boundary, the bulk free energy Fj, is minimized 
using standard techniques of the calculus of variations. 
In the case of strong anchoring it is also sufficient to 
minimize the bulk free energy but subject to n taking 
prescribed fixed values at the boundary. In the case of 
weak anchoring the total free energy F, which includes 
the surface free energy F S) has to be minimized. The 
most commonly used expression for the surface free en- 
ergy is of the form proposed by Rapini and Papoular by 



F s [n(r) 



1 



d 2 r w(r)(h ■ u) 



(3) 



The integral runs over the boundary and w = w(r) is the 
corresponding anchoring strength that characterizes the 
surface. The local unit vector perpendicular to the sur- 
face is denoted as n. For negative w, this contribution fa- 
vors an orientation of the molecules perpendicular to the 
surface, while positive w favor degenerate planar orienta- 
tions at the surface. The absolute value of the anchoring 
strength is commonly of the order 10~ 6 to 10 -2 N/m. 



B. The model 

Here we consider a nematic liquid crystal confined be- 
tween a patterned substrate at z = and a flat substrate 
at z = D where the z axis is normal to the flat substrate. 
As Figs. [1] (a), (c), and (d) illustrate, the lower substrate 
is characterized by geometrical and/or chemical patterns 
of periodicity p along the x axis. Moreover, the system 
is translationally invariant in the y direction. Within 
the one-constant approximation [Eq. ([2])] the total free 
energy functional is given by 

F[0(x,z)',e D ,D] = 

p D 

^ j ' dx J dz[V6(x,z)¥+F s [6{x,z {x))],{l) 

z a {x) 

where L is the extension of the system in y direction, 
z (x) is the surface profile of the patterned substrate, and 
n = (sin0, 0, cos#). At the upper surface strong anchor- 
ing 6(x, z = D) = 9d is imposed. Twist is not considered, 
i.e., K22 — 0. Of course, the analysis can be straightfor- 
wardly extended to the case of different splay and bend 
constants Ku and ^33, respectively, but this aspect of 
the problem is important only if the analysis is supposed 
to yield quantitative results for a specific nematic liquid 
crystal. The surface contribution F s [9(x, zq(x))] includes 
the anchoring energy for which the Rapini-Papoular form 
[Eq. ©] is adopted: 
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FIG. 1: (a) The system under consideration consists of a 
nematic liquid crystal confined between two substrates at a 
mean distance D. The upper flat substrate induces strong 
anchoring, i.e., 9(x, z = D) = 6d, while the lower substrate is 
characterized by a surface pattern of period p in x direction. 
The system is translational invariant in y direction perpendic- 
ular to the plane of the figure. In (a) the lower sinusoidally 
grating surface (groove depth A) is endowed with an alter- 
nating stripe pattern of locally homeotropic anchoring (white 
bars) and homogeneous planar anchoring (black bars). The 
period of the chemical pattern is half the period of the sur- 
face grating p. In (c) and (d) a pure geometrically struc- 
tured lower substrate and pure chemically patterned lower 
substrate, respectively, are shown. The anchoring direction 
at the substrates is schematically represented by black rods. 
Quantitatively reliable predictions of the phase behavior of a 
nematic liquid crystal confined between two substrates at an 
arbitrary mean distance D in (a) can be achieved if the effec- 
tive free energy function [Eqs. © - ©] i s analyzed. To this 
end the free energy functional of the nematic liquid crystal 
confined between two substrates at a single and rather small 
mean distance D* in (b) has to be minimized as is discussed 
in the main text. 



F s [9(x,z (x))} = 
— r J .... ( . ) (-Bin(go(x)K(x) + cos(ffl (x))) 2 
2 Jo V / l+» 

where 9q(x) = 9{ (x)). Equations @ and (J5J 

together with the boundary condition for the nematic di- 
rector at the upper surface completely specify the free en- 
ergy functional for the system under consideration. The 
Euler-Langrange equations resulting from the stationary 
conditions of the total free energy with respect to the 
nematic director field can be solved numerically on a 
sufficiently fine two-dimensional grid using an iterative 
method. In order to obtain both stable and metastable 
configurations, different types of initial configurations are 
used. However, due to the pattern of the lower surface 
the determination of the director field and the phase dia- 
gram turns out to be a challenging numerical problem in 
particular in the case of large cell widths D. Moreover, 



the energy barrier between two metastable states cannot 
be determined this way. 



C. Effective free energy function 

Here we map the free energy functional 
F[8(x, z); Ojj, D] [Eq. Q] of a nematic liquid crys- 
tal cell with arbitrary width D and arbitrary anchoring 
angle 9d at the upper surface (see Fig. Q] (a)) onto the 
effective free energy function 

fW(0 o ,o d ,d) = S (e D - e o y+F^(e ),(6) 

where the average surface director orientation 6> 62] at 
the lower patterned surface is given by 

D* 8 

0o(#d.,-D*) = Ob, — — F(0D t ,D*)\ m i n . (7) 



KLp 86 D , 



The effective surface free energy function Fi e charac- 
terizing the anchoring energy at the patterned surface 
can be written as 

F^f\e (e D „D*),D*) = 

F{9d„ , D*)\ m i n — — — — (do, — §q(9d, j D *)") . 



(8) 



In order to calculate 9q and Fs e ^^ explicitly, we first 
determine numerically the minimum of the free energy 
F(9d,, D*)\min of the nematic liquid crystal cell for a 
single and rather small value D = and arbitrary an- 
choring angle 9d, at the upper surface (see Fig. Q] (b)). 
Thereafter the phase behavior, energy barriers between 
metastable states, and effective anchoring angles for the 
system of interest (Fig. [T] (a)) can be obtained for ar- 
bitrary values of D and 9d from F^^(6q, 9d, D) as a 
function of the single variable #o- Such a calculation is 
considerably less challenging than minimizing the origi- 
nal free energy functional F[9(x, z);9d, D] with respect 
to 9(x, z) on a two-dimensional (x, z) grid. However, the 
effective free energy method is applicable only if Eq. ([7]) 
can be inverted in order to obtain On, (9o) which is needed 
as input into Eq. ([8]) . The condition for this inversion fol- 
lows from Eq. ([7]): 



£>. d 2 
KLp 89% 

Aeff) 



■F(e D „D.)\, 



> 0. 



(9) 



Moreover, Fs (0o(#d» ,D*),D*) is practically indepen- 
dent of the cell width provided > p implying that 
the interfacial region above the lower substrate does not 
extend to the upper substrate. 

Before studying the nematic liquid crystal in contact 
with the patterned substrates shown in Fig. [T] it is in- 
structive to analyze first the nematic liquid crystal con- 
fined between two homogeneous fiat substrates at a dis- 
tance D. The upper surface induces strong anchoring, 
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i.e., 9(z = D) = Oo- The free energy functional defined 
in Eq. (g]) follows as 

F{9(z);0 D ,D] = ^f jdz {^j + F S (9 Q ).(10) 
o 

The solution of the Euler-Langrange equation d^9(z) = 
subject to the boundary condition at the upper surface 
9(z = D) = Ojj interpolates linearly between the top and 
bottom surfaces: 



s{z) = e D --{D- z ){e D -e ), (11) 

where 9{z = 0) = 9q follows from the boundary condition 
at the lower surface. With this solution of the Euler- 
Lagrange equation the minimized free energy function 
reads 

F(0D« , A.)|min = -^j^ {®D, ~ 6of + F s (6q) . 

(12) 

It follows directly from Eqs. © - jH) that 6 = $o, 
F s (e//) = F s , and 

F^f\6 ;6 D ,D) = ^ (9 D - 9 ) 2 + F s (9 ). (13) 

Hence in the case of homogeneous confining substrates 
the effective free energy function [Eq. (TT3|) ] agrees exactly 
with the minimized free energy function of the original 
system [Eq. (fT2j) ]. 

III. APPLICATIONS 

A. Geometrically and chemically patterned 
substrates 

In the previous section we have shown that we can 
describe a nematic liquid crystal confined between two 
homogeneous substrates by the effective free energy func- 
tion [Eq. ([B]) ], In this subsection we apply this approach 
to the particular case of a nematic liquid crystal confined 
between a chemically patterned sinusoidal surface and 
a flat substrate with strong homeotropic anchoring (see 
Fig. [1] (a)). The surface profile of the grating surface is 
given by zq(x) = Asm (qx), where A is the groove depth 
and p = 2n/q is the period. As Figure [1] (a) illustrates, 
the surface exhibits a pattern consisting of alternating 
stripes with locally homeotropic and homogeneous pla- 
nar anchoring. The projection of the widths of the stripes 
onto the x axis is p/A and the anchoring strength is spec- 
ified by a periodic step function: w(x) — —wh and wp 
for values of x on the homeotropic and planar stripes, 
respectively. Figure (a) displays F(9d, 7 D^ = p)\ m in 
(dashed line) and i 71 ] 6 ^ (#o, -D* = p) (solid line) for 



pw H /K = 1, pwp/K = 2.5, and A/p = 0.09. The 
shapes of F(9p) itl D :t = p)\ m i n as a function of 9p,, and 
Fs £ ^\9~o, = p) as a function of 9q are rather similar 
because 9o ~ #d» (see Eq. (JSj) ) for this set of model pa- 
rameters. Figure [2] (b) displays the phase diagram plot- 
ted as a function of the anchoring angle 9jj at the upper 
substrate and the mean separation of the substrates D. 
The calculations demonstrate the existence of two (sta- 
ble or metastable) nematic director configurations: the 
homeotropic (H) phase, in which the director field is al- 
most uniform and parallel to the anchoring direction im- 
posed at the upper surface, i.e., njy = (sin#D, 0, cos#d), 
and the hybrid aligned nematic (HAN) phase, in which 
the director field varies from n# at the upper surface 
to nearly planar orientation through the cell. Note that 
there are two HAN textures: HAN + and HAN_ corre- 
sponding to positive and negative average surface angles 
at the lower surface (see also Fig. [7] below). For small 
anchoring angles 9d the HAN phases are stable provided 
the cell width is larger than D coex (more precisely, the 
HAN+ texture is stable for 9 D > while the HAN_ tex- 
ture is stable for 9d < and they coexist at 9d = 0). 
For smaller distances between the substrates D < D coex 
the HAN phases are no longer stable because distortions 
of the director field are too costly in the presence of the 
dominating strong anchoring at the upper surface. The 
comparison of the phase boundary of thermal equilib- 
rium as obtained from the effective free energy method 
[Eqs. © - ©, and solid line in Fig. [U (b)] and the di- 
rect minimization of the underlying free energy functional 
[Eqs. (0| and ©, and diamonds in Fig. [2] (b)] demon- 
strate the reliability of the effective free energy method. 

We note that the phase transition between the H and 
HAN textures is first order despite the fact that the effec- 
tive surface free energy favors monostable planar anchor- 
ing, i.e., Fg e ^^ exhibits only a minimum at 9q = ir/2 in 
the interval 9q € [0,7r/2]. A first order phase transition 
in a nematic liquid crystal device with a monostable an- 
choring condition on a homogeneous lower substrate has 
been predicted for the special case 9p, — in Refs. [4^.l48j 
using the empirical expression 

F s (9 ) = ^sin 2 (20 o ) + Wl sin 2 (0 o ) (14) 

as input into Eq. (I13| . For wi < 2wq this surface free 
energy has two minima at #o = and 9 = n/2 in the 
interval #o € [0, 7r/2], and as such is bistable, while for 
w\ > 2wq, only the minimum 9q ~ exists, i.e., the 
surface is monostable. To study the stability limit of the 
H phase we expand F < - eff ^(9 ) in Eq. §T5§ around 9 = 
up to sixth order: 

+ ^(0)9^ + ^(0)91 (15) 
The H phase corresponds to a local minimum of 
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FIG. 2: (a) The minimized free energy F(8o,,D* = p)\ min 
(dashed line) and the effective surface free energy 



FIG. 3: (a) The minimized free energy F(8d,,D* = p)\ min 
(dashed line) and the effective surface free energy 
Fr JJ '(e a , D, = p) (solid line) of a nematic liquid crystal con- p Wf) Dt = p ) ( solid line ) of a nema tic liquid crystal con- 



fined between a flat surface with strong anchoring and a 
chemically patterned sinusoidal surface with groove depth 
A/p = 0.09 (see Fig. Q] (b)). L is the extension of the cell 
in the invariant y direction, K is the isotropic elastic con- 
stant, and the anchoring strength on the homeotropic stripes 
(white bars in Fig.[T](b)) and planar anchoring stripes (black 
bars in Fig.[l](b)) are pwn/K = 1 and pwp/K = 2.5, respec- 
tively, (b) Phase diagram of the same system as a function 
of the anchoring angle at the upper flat substrate 8d and the 
cell width D (see Fig. Q](a)). The solid line denotes first or- 
der phase transitions between a homeotropic (H) and hybrid 
aligned nematic (HAN+ and HAN_) phases. At 9d = and 
D/p ~ 5.8 there is a triple point where the HAN+, HAN_, 
and H states coexist. The solid circle marks the critical point 
at Dcr/p ~ 3 and #^ r ' ~ ±15°. The limits of metastabil- 
ity of the HAN+ (1) and the H (2) state are denoted by the 
dot-dashed lines. The limit of metastability of the H state for 
9d < and the HAN_ state are not shown for clearness. The 
lines and the solid circle follow from analyzing the effective 
free energy function [Eqs. ((6]) - (JSj>] while the diamonds repre- 
sent the phase boundary of thermal equilibrium as obtained 
from a direct minimization of the underlying free energy func- 
tional [Eqs. (0J and ©]. 



fined between a flat surface with strong anchoring and a 
chemically patterned sinusoidal surface with groove depth 
A/p — 0.09 (see Fig. [1] (b)). The anchoring strength on 
the homeotropic stripes and planar anchoring stripes are 
PWh/K — 2.5 and pwp/K — 6, respectively, (b) Phase 
diagram of the same system as a function of the anchoring 
angle at the upper flat substrate 6d and the cell width D (see 
Fig.[T](a)). The solid line denotes the first order phase transi- 
tion between a homogeneous (H) and hybrid aligned nematic 
(HAN+) phase. The solid circle marks the critical point at 
w 0.4 and 6»^ r) w 29°. The limits of metastability of 
the HAN+ (1) and the H (2) state are denoted by the dot- 
dashed lines. The lines follow from analyzing the effective 
free energy function [Eqs. - (JSj) ] while the diamonds and 
the solid circle represent the phase boundary and a critical 
point as obtained from a direct minimization of the underly- 
ing free energy functional [Eqs. Q and (J5}]. For D < p the 
phase transition as well as the limits of metastability cannot 
be determined using the effective free energy function because 
Fs (0o) is not known in the region close to its maximum. 
For clearness, only the phase diagram for positive Qd is shown 
(c.f. Fig.rj. 



F^f \e Q ) in Eq. p]) if D < D^. where 

KLp 



(16) 



A standard bifurcation analysis reveals that the transi- 
tion from the H phase to the HAN phase can be either 
first order or continuous. The transition is continuous if 



Fj 4) (0) > 0, first order if F s (4) (0) < 0, and F s (4) (0) = 
corresponds to a tricritical point. In the case of a first 
order phase transition the phase boundary of thermal 
equilibrium is given by 



10 (Fj 4) (0) 



Dtr 16Fj 6) (0) sW 



(17) 
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We emphasize that the order of the phase transition de- 
pends only on the surface free energy F s (9q) close to 
do = for 0£> = 0. Therefore it is possible to have a first 
order phase transition even with a monostable surface 
characterized by a monotonic surface free energy such as 
the one shown in Fig. [2] (a) for 9o £ [0, 7r/2] or the empiri- 
cal equation (114|) with w\ > 2wo as well the more general 
expression 

oo 

F s (0 ) = ^ [a„ cos(2n6» ) + 6„ sin(2n6» )] (18) 

with appropriate parameters a n and b n . 

First order phase transitions between the H and HAN 
texture are of particular interest for bistable liquid crys- 
tal displays. In a bistable liquid crystal display the 
two molecular configurations corresponding to light and 
dark states are locally stable in the therm ody namic space 
when the applied voltage is removed [U|5J|- Therefore, 
power is needed only to switch from one stable state to 
another, in contrast to monostable liquid crystal displays 
which require power to switch and to maintain the light 
and the dark states. 

We now turn our attention to the case that it is 
not possible to evaluate 0d» (#o,£>*) from #o(#d„ , -D*) 
[Eq. ([7])] because the condition for this inversion is not 
satisfied [Eq. ©]. To this end we have chosen the pa- 
rameters pwn/K = 2.5, pwp/K — 6, and A/p = 0.09 
for the system shown in Fig. [1] (a) . Figure [3J (a) display 

F(0 D . , D* = P)\min (dashed line) and F^ ff) (9 , D* = p) 
(solid line) while the corresponding phase diagram is 
shown in Fig. [3J (b). As is apparent from the solid line in 
Fig.[3J(a) it is not possible to determine the effective sur- 
face free energy function for the all values of 9a because 
the upper flat substrate at D* is too far away from the 
lower patterned substrate in order to induce all possible 
average anchoring orientations 9q- I n other words, the 
anchoring energy at the patterned substrate is too large 
to be balanced by the elastic energy for the chosen mean 
distance Z?» between the substrates (see Fig.[T](b)). Nev- 
ertheless Fig. [3J (b) demonstrates that even this partial 
information about the effective surface free energy func- 
tion can be used to calculate the phase diagram for cell 
widths sufficiently larger than the width at the critical 
point D cr . 



B. Purely geometrically structured substrates 

In the last subsection we have shown that with a suit- 
able chemical and geometrical surface morphology on one 
of the interior surfaces of a liquid crystal cell, two stable 
nematic director configurations can be supported. The 
zcnithally bistable nematic devices that have been stud- 
ied recently [46|, H3, IH, |4^| consist of a nematic liquid 
crystal confined between a chemically homogeneous grat- 
ing surface (see Fig.Q](c)) and a flat substrate with strong 



homeotropic anchoring. The profile of the asymmetric 
surface grating is given by 

zq[x) — As\n.(qx + hsm(qx)) , (19) 

where A is the groove depth, p — 2n/q the period, and 
h is the "blazing" parameter describing the asymmetry 
of the surface profile. Such a grating surface has been 
studied by Brown et al. [l5| who found a first order 
transition between the HAN state, characterized by a low 
pretilt angle (9o ~ an( i the H state, characterized 

by a high pretilt angle (9q ~ 0). Strictly speaking, the 
H state does not correspond to the homeotropic texture 
(see Fig. [S](b) below), but we keep the same notation as 
in the previous section for consistency. Here we study 
phase transitions of a nematic liquid crystal in contact 
with the blazed surface in a more detail using the effective 
free energy method discussed in Sec. II C. 

In Fig. d] (a) the minimized free energy F(9r>, , -D* = 
p)\min (dashed line) and the calculated effective surface 
free energy Fs e ^\§o, £>* = p) (solid line) are shown for 
the anchoring strength pwu / K = 2 on the grating sur- 
face, the groove depth A/p = 0.27, and the blazing pa- 
rameter h = 0.2. The effective surface free energy is 
asymmetric with respect to 9q = because of the asym- 
metry of the grating surface. As a consequence also the 
phase diagram, plotted as a function of the anchoring 
angle on the upper surface 9d and the distance D/p is 
asymmetric (Fig. [4] (b)). However, the topology of the 
phase diagram is the same as in the case of a symmetric 
substrate (see Figs. [U (b), [3J (b), and for a more general 
discussion Sec. Ill C below). 

We now concentrate on the most interesting (from a 
practical point of view) case of strong homeotropic an- 
choring (8o = 0) at the upper homogeneous surface. Fig- 
ure [5] (a) displays the phase diagram for a few values of 
the blazing parameter h and a fixed value of the local 
homeotropic anchoring strength on the grating surface 
pwn/K = 2. For a fixed cell width D/p, asymmetry 
(h ^ 0) leads to a decrease of the groove depth A/p at 
which there is a first order transition between the HAN 
and the H phases, as compared to the nematic liquid 
crystal cell with the symmetric surface grating [h = 0). 
Upon increasing the groove depth A/p the transition line 
ends at a critical point (not shown in the figure), while it 
diverges as A — ► A Q . The groove depth A corresponds 
to an anchoring (or surface) transition between low tilt 
and high tilt surface states which are the homeotropic 
and planar effective anchoring states, respectively, in the 
case h = 0. 

The effective free energy method allows one to calcu- 
late an energy barrier E between the two bistable states, 
which is not feasible by the direct numerical minimiza- 
tion of the free energy functional. The results of our 
calculations are shown in Fig. O (b). The asymmetry of 
the surface grating leads to a decrease of the energy bar- 
rier. With increasing groove depth A/p the energy bar- 
rier decreases and eventually vanishes upon approaching 
the critical point. The diamonds in Fig. [U(b) denote the 
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FIG. 4: (a) The minimized free energy F(8d,,D* = p)\min 
(dashed line) and the effective surface free energy 

Fa {0o, D t = p) (solid line) of a nematic liquid crystal con- 
fined between a flat surface with strong anchoring and a chem- 
ically uniform blazed surface (see Fig.[T](c) and Eq. (|19p ). The 
locally homeotropic anchoring strength on the blazed surface 
is pwu/K — 2, the groove depth is A/p = 0.27 and the blaz- 
ing parameter is h = 0.2. L is the extension of the cell in 
the invariant y direction and K is the isotropic elastic con- 
stant. The effective surface free energy Fj e ^' (as well as 
F(8d«,D* = p)\min) is periodic with the period 2-n but it is 
asymmetric with respect to #o = 0. (b) Phase diagram of the 
same system as a function of the anchoring angle at the upper 
flat substrate 9d and the cell width D with the same line code 
as in Figs.[2]and[3] The triple point (where the HAN+, HAN_ 
and H phases coexist) is at 6d ~ 12° and D/p » 8, and the 
critical points (solid circle) are at D cr /p « 1.3, 9^ r) « 38° 



and D C r/p ~ 1-7, 9 



(cr) 



-24°. For clearness, only the limits 



of metastability of the HAN_ (1) phase and the H (2) phase 
(for negative 9d) are shown. 



energy barriers for the abovementioned anchoring transi- 
tions of a nematic liquid crystal in contact with a single 
grating surface. 

The energy barrier between two bistable states is an 
important quantity for the design of a zenithally bistable 
nematic device. Too small energy barrier, as compared 
to ksT, would cause spontaneous switching between the 
two states because of thermal fluctuations, while en- 
larging the energy barrier leads to an increase of the 
power consumption. Using the calculated values of E 
(see Fig. [5] (b)) one can estimate the energy barrier in a 
real nematic liquid crystal cell. For instance, for a cell 
of area 1 /xm x 1 /jm and of width D = 5 /xm, and taking 
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FIG. 5: (a) The phase diagram of a nematic liquid crys- 
tal confined between the blazed surface (see Fig. [T] (c) and 
Eq. (|19p ) with local homeotropic anchoring and a flat surface 
with strong homeotropic anchoring (9d = 0) as a function of 
the groove depth A/p and the cell width D/p. The anchor- 
ing strength on the blazed surface is pwn/K — 2. The lines 
correspond to different values of h and denote first order tran- 
sitions between homeotropic (H) and hybrid aligned (HAN_) 
phases. For small A/p, the lines extend to D = oo corre- 
sponding to a first order anchoring transition between planar 
and homeotropic phases. Upon increasing A/p the first order 
transition lines end at critical points which are not shown in 
the figure, (b) The energy barrier at the first order transitions 
with the same line code as in (a). The lines have been ob- 
tained from the total effective free energy (see Eq. (0), while 
the diamonds correspond to the energy barriers between the 
planar and homeotropic effective anchoring which follow from 
considering only the surface contribution Fs e!! \ 



the typical values K = 5 x 1CT 12 N and wh = 1CT 5 N/m, 
one obtains E w 37 k B T for h = 0.2 and A w 0.026 //m, 
which seems to be an acceptable value. 

Another important quantity in zenithally bistable ne- 
matic devices is the average director orientation at the 
grating surface in the two degenerate states. The average 
surface director in the HAN (0q HAN ' > ) and H ) states 
is shown in Fig. [6] for the same model parameters as in 
Fig. [5] and for the values of A/p and D/p corresponding 
to the coexistence line. The asymmetry of the surface 
grating leads to a decrease of 0q HAN ' ) and an increase of 
6>q . For a fixed value of h, the difference between the 
two angles decreases with increasing the groove depth 
A/p and finally vanishes upon approaching the critical 
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FIG. 6: The average surface director (see Eq. (0)) for the 
hybrid aligned (HAN_) phase in (a) and the homeotropic 
(H) phase in (b) on the lines of the first order transitions 
(see Fig. for a nematic liquid crystal confined between 
the blazed surface (see Fig. [1] (c) and Eq. (|19[l ) with lo- 
cal homeotropic anchoring and a flat surface with strong 
homeotropic anchoring (9d = 0). The diamonds in (a) and 
(b) correspond to the two minima of the effective surface en- 
ergy function Fs (9o). The model parameters and the line 
code are the same as in Fig. [5] 



point (not shown in the figure). 

Hence the asymmetry of the surface grating leads to 
a decrease of the groove depth at which the Instabil- 
ity is observed, which improves optical properties [63| . 
and to a decrease of the energy barrier, which lowers the 
power consumption of a zenithally bistable nematic de- 
vice. On the other hand, also the difference between the 
two bistable states decreases which impairs optical prop- 
erties of such a device. 



C. Phase diagrams for a model surface free energy 



In section IIII Al we have discussed phase diagrams in 
the (9d,D/p) plane which can be described in terms of 
the surface free energy given by Eq. (|14[) . However, it 
is instructive to consider a more general situation that 
the surface free energy follows from a truncation of the 
Fourier expansion given in Eq. {T5J . To be able to study 
both symmetric and asymmetric surfaces we assume a 



natural generalization of Eq. (|14[) , namely 

2 

F„(9 ) = K cos(2n6» ) + b n sin(2n6> )] , (20) 

71=0 

which reduces to Eq. (fT4")) in the case of a symmetric 
surface characterized by b n = 0. The angle #( mm ) that 
minimizes F^^(0q] #d, D) (see Eq. (fT5|) ) is a function 
of Op and D. The derivative 
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(21) 



is the susceptibility of the system that diverges at the 
critical thickness 



D cr = -KLp/ mm F?(9 ), 

So 



(22) 



and remains finite and positive for D < D cr . From 

Eq. ([22]) the conditions for the critical angle 9q follow 
as 

Fs (3)^(cr) ) = and F W(0(cr)) >o> (33) 

implying that D cr and 9^ T ^ depend only on the form of 
F s (6 ). 

The extremes of F s given by Eq. (|20|) can be found 
easily only in the case of symmetric or antisymmetric 
(a n — 0) surface and the same concerns the position of 
critical point. In this work, however, we are interested 
rather in possible topologies of the phase diagram in the 
{9 D ,D/p) plane, which result from Eqs. (JT3J) and pT]) . 
and not in the exact location of critical points or tran- 
sition lines. To draw schematic phase diagrams we con- 
sider F s as a function defined on the unit circle z = e 2l6 ° . 
Depending on the parameters a n and b n , F s (9q) has ei- 
ther one minimum and one maximum or two minima and 
two maxima. This conclusion applies also to the function 
F"(0o), thus, there can be either one or two critical points 
in the phase diagram (see Eqs. ([22]) and (|23|) h In the 
limit of large D, there is always a first order phase transi- 
tion between two non-uniform textures corresponding to 
the opposite orientations of the director at z = Q. This 
is because F^ e ^^ is not a periodic function of 9q at fixed 
9d- To find 9o at the transition we expand F s around 
its deepest minimum (denoted 9 m ), which leads to the 
approximate free energy: 



F^f\9 D ,D)\ r 



F, 



+ 



KLp(9 D - 9 m f 
2(D + b) 



(24) 



where b = KLp/ F'J (9 m ) is the extrapolation length. 
Since 9 m and 9 m ± t are equivalent minima of F s , and 
both 9 m and 9]j are allowed to vary in the interval 
[— 7r/2, 7r/2], the transition occurs at 6^ r ' = 9 m + n/2 if 



-7r/ 2 < #m < or at 9^ r> = 9 m - tt/2 if < 9 m < tt/2. 
With the above information we can now draw schemat- 
ically the phase diagram (see Fig. [7]) . Since the vertical 
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FIG. 7: Schematic presentation of possible types of the phase 
diagram in the (8d,D/p) plane for the surface free energy 
given by Eq. (|20|l . To and T± denote different (usually non- 
uniform) textures, the thick lines correspond to first order 
transitions, and black circles mark critical points. When F a 
(defined on the unit circle z = e 2l8 °) has one minimum and 
one maximum the phase diagram can be either of type (a) or 
(b). When F s has two minima and two maxima the phase 
diagram is of type (b) unless the minima are of equal depth, 
in which case it is of type (c). Note that the lines 8d = in/2 
are identified with each other and the phase diagram can be 
considered as being on a cylindrical surface. 



lines at 9d = ±7r/2 are to be identified with each other 
the phase diagram can be considered on a cylindrical sur- 
face. Away from the transition lines there is a smooth 
evolution from one texture to another. This means that 
at fixed D it is possible to transform smoothly the T + 
texture into the T_ texture, i.e., without crossing the 
transition line, even for D > D cr . We note that in some 
range of parameters, the phase diagram for F s with one 
minimum is topologically indistinguishable from that for 
F s with two minima (Fig. (b)); in both cases there are 
two critical points and a triple point. If F s has two equal 
minima (e.g., at 6q = and 9q = ±w/2 in the case of 
symmetric surface) the triple point disappears and the 
first order transition lines extend to D — > oo, as shown 
in Fig. [7] (c). 



IV. SUMMARY 

We have studied the phase behavior of a nematic liquid 
crystal confined between a flat and a patterned substrate 
(Fig. 0) using the Frank-Oseen model [Eq. (]3|)] and the 



Rapini-Papoular surface free energy [Eq. ([5])] . An expres- 
sion for the effective free energy function of the system 
[Eq. ©] was derived by determining an effective surface 
free energy characterizing the anchoring energy at the 
patterned surface [Eq. ©]. Using the effective free en- 
ergy function, we have determined the phase behavior of 
the nematic liquid crystal confined between a fiat surface 
with strong anchoring and a chemically patterned sinu- 
soidal surface (Fig. Q] (a)), finding first order transitions 
between a homeotropic texture (H) and hybrid aligned 
nematic (HAN) textures (Figs. (b) and 0(b)). It is 
possible to have a first order phase transition even with 
a monostable surface characterized by a monotonic sur- 
face free energy function (Fig. (a), 9 £ [0,7r/2]). In 
addition we have performed direct minimizations of the 
original free energy functional [Eqs. (0| and (0)] on a two- 
dimensional grid and found remarkably good agreement 
with the phase boundaries resulting from the effective 
energy function analysis (Figs. (b) and 0(b)). Hence 
quantitatively reliable predictions of the phase behavior 
can be achieved using the effective free energy method. 

Using this method, we have also studied the phase be- 
havior (Fig. [31(b)) of a nematic liquid crystal confined be- 
tween a chemically uniform, asymmetrically grooved sub- 
strate (Fig. 0(c) and Eq. (|19[) ) with locally homeotropic 
anchoring and a flat substrate with strong homeotropic 
anchoring, which is a typical setup for a zenithally 
bistable nematic device [ifl 53, El, S§| • The asymmetry 
of the grating substrate leads to a decrease of the groove 
depth at which a first order transition between the H and 
HAN phases occurs (Fig. [5] (a)). Moreover, we have de- 
termined the energy barrier between the two coexisting 
states (Fig. (b)). Our calculations show that the en- 
ergy barrier decreases with increasing the asymmetry of 
the grating surface but it is well above ksT for a typical 
nematic liquid crystal cell. In addition, the average di- 
rector orientation at the grating surface in two bistable 
states has been calculated (Fig. [5]). The difference be- 
tween the two bistable states vanishes with increasing 
substrate asymmetry, which has a negative effect on the 
optical properties of a zenithally bistable nematic device. 

We have also generalized the model of the effective sur- 
face free energy considered by Parry- Jones et al. [H, |48| 
to the case of asymmetric structured substrates and ob- 
tained three possible types of the phase diagram in the 
plane spanned by the orientation of the director at the 
homogeneous surface and the thickness of the nematic 
cell. The asymmetry of the substrate causes only a shift 
of transition lines and critical points, compared to the 
symmetric case, but does not change the topology of the 
phase diagram. Finally, we have verified that this model 
allows one to reproduce qualitatively the phase diagram 
of a nematic liquid crystal confined between a homoge- 
neous planar substrate and an asymmetrically grooved 
surface (Fig. (b) and Fig. 0(b)). 
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